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t^> ', Abstract 

The fireball created in an ultrarelativistic heavy ion collision is the en- 
vironment in which all processes providing clues about the possible for- 
mation of the quark-gluon plasma (QGP) happen. It is therefore crucial 
' to understand the dynamics of this hot and dense system. We set up 

0^ . a model in which the fireball evolution is reconstructed between two 

, stages, the freeze-out, which is accessible by hadronic observables, and 

£NJ ' the initial collision for which the overlap geometry can be calculated. 

Using the equation of state (EoS) provided by a quasiparticle model 
. of the QGP, we are able to calculate thermodynamical properties in 

volume slices of constant proper time and determine the volume expan- 
O ■ sion self-consistently. The resulting evolution model can then be tested 

q^. against other observables, such as dilepton yields. 

' 

1 Introduction 

•i-H . 

^ ■ One of the most striking predictions of QCD at high temperatures is the for- 

mation of the quark-gluon plasma (QGP), a deconfined phase with quarks and 
gluons as degrees of freedom. Unfortunately, finding direct evidence for this 
phase or studying its properties is not an easy task. Ideally, one would like to 
study the behaviour of an observable sensitive to the formation of the QGP, 
such as the J/\l/ system or dilepton radiation. In reality, however, no static 
QGP system can be prepared in a heavy ion collision and all processes inside 
the resulting fireball are necessarily convoluted with its evolution dynamics. 
As this leaves two a priori unknown inputs into any model calculation, the 
modified reaction dynamics inside the system and the fireball evolution itself, 
a model which is able to describe just one (or a few) observables may not be 
sufficient to disentangle these two pieces and extract unambiguous information. 



1 Work supported in part by BMBF and GSI. 



The proper way to proceed seems therefore as follows: From the theoretical 
side, and especially from lattice calculations (see e.g. |J), much is known 
about thermodynamical properties of the QGP in a static situation. Although 
there is no direct possibility to compare to experiment, this knowledge can be 
used under certain circumstances as an input for dynamical fireball models, 
which constitute the next important building block for modelling any specific 
process inside such a system. Here, experimental input is available in the 
form of bulk hadronic observables such as directed flow and total multiplicities 
which characterize the evolution endpoint. The third and final step is then 
the calculation of a given process within the arena prepared by the fireball 
evolution. In this way, fireball evolution and the process in question are tied 
to different observables and are readily disentangled. 

In the following, we try to elaborate on the second building block, the 
dynamical fireball evolution in a simple framework and calculate as an appli- 
cation dilepton radiation from this fireball. This is more completely worked 
out in 0. 

2 Model framework 

We do not aim at the detailed description of the heavy-ion collision or its 
subsequent expansion on an event by event basis or by a hydrodynamical 
simulation. Instead we use a model of an expanding fireball which enables us 
to rapidly test different scenarios with different parameter sets in a systematic 
way, so as to gain insight into the time evolution of the strongly interacting 
system. 

We assume that the physics of the fireball is the same inside each volume 
of given proper time r, thus averaging over spatial inhomogenities in density 
and temperature. The volume itself is taken to be an expanding cylinder, in 
which the volume elements move away from the center in order to generate the 
observed flow. There is no global Lorentz frame in which thermodynamics can 
be applied. As the fireball expands, volume elements away from the center are 
moving with large velocities and are subject to time dilatation when seen in 
the center of mass frame of the collision. In this frame, the fireball expands, 
at a given time, much more rapidly in the center than at the edges and does 
not resemble a cylinder any more. We assume a linear increase in rapidity 
when going from central volume elements to the fireball edge along the beam 
(z)-axis and the transverse axis. As the velocities along the z-axis are typically 
large (up to c) as compared to transverse motion (up to 0.55 c) for SPS and 
RHIC conditions, we make the simplifying assumption that the proper time 



is in a one-to-one correspondence to the ^-position of a given volume element, 
thus neglecting the time dilatation caused by transverse motion. The whole 
system is assumed to be in local thermal (though not necessarily chemical) 
equilibrium at all times. 

Given this overall framework, the volume expansion of the fireball is gov- 
erned by the longitudinal growth speed v z and the transverse expansion speed 
vj_ at a given proper time. These quantities can be determined at the freeze- 
out point and correspond to the observed amount of flow. However, flow is 
measured in the lab frame and needs to be translated into the growth of proper 
time volume. 

Note that it is important to keep track of the velocity v z at the fireball 
edge in this setup. If this quantity is simply assumed to be c, the volume of 
given proper time does not grow like rc as one might naively expect (the cm. 
frame volume does so, however), but is infinite right from the beginning. 

We use a detailed analysis of the freeze-out conditions for central Pb-Pb 
collisions at 160 AGeV P] to fix the endpoint of the evolution. The initial state 
is constrained using the overlap geometry of the colliding nuclei. The expansion 
between initial and freeze-out stages is then required to be in accordance with 
the EoS as determined from the quasiparticle model described in ||. 

The volume expansion is parametrized by the following set of equations: 

«l(t) = f dr'c ± ^)- R(t) = R + f f dr'dr"c ± ^Q (1) 
Jo e (r) Jo Jo e(r ") 



VzM = < + f dt'cj^)- z(t) = z + vl-t+ f f dUdfc^Q-. (2) 
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Here, the acceleration was assumed to be proportional to the ratio of pres- 
sure p over energy density e with a proportionality constant c. The free pa- 
rameters c_i_,c 2 , freeze-out proper time 77 and freeze-out cm. time tf can be 
fitted by requiring agreement with initial conditions Rq ~ 4.5 fm and v ± = 
(overlap geometry) and final conditions Rf ~ 8.55 fm, v±_ = 0.5c, Tf = 100 
MeV, v z = 0.9c (results from 0). Assuming entropy conservation with an 
entropy per baryon of 26 for SPS conditions at 160 AGeV, the entropy density 
s at a given proper time can then be obtained by dividing the total entropy 
So by the volume, s = S /V(t). With the help of the EoS, the temperature 
T(s), pressure p(s) and energy density e(s) can then be calculated. These are 
inserted into eqs.(l-2) in order to yield a self-consistent solution. Chemical 



potentials for hadrons are introduced in order to agree with the experimen- 
tally observed particle abundancies. For the above conditions, a rise in the 
pion chemical potential fj, n up to 123 MeV towards freeze-out is necessary. 



3 Results 

The resulting volume expansion and temperature profile for SPS conditions at 
two different energies is shown in Fig. |l| 
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Figure 1: Left panel: Volume of constant proper time at SPS 160 and 40 AGeV 
for central Pb-Pb collisions. Right panel: Corresponding temperature profile 
with entropy per baryon s/pb = 26 (160 AGeV) and 13 (40 AGeV). 



The primary uncertainty is attached to the initial temperature, which ap- 
pears rather large (~300 MeV) for 160 AGeV collisions. Part of the difference 
in this value as compared to other approaches comes from the use of a more 
realistic equation of state, which even for the initial temparatures differs con- 
siderably from the EoS of an ideal gas. Moreover, it depends on the initial 
longitudinal expansion velocity. After a few fm/c evolution time, however, 
the system is not very sensitive to details of the initial conditions any more. 
Remarkably, even the 40 AGeV scenario shows some evolution above the crit- 
ical temperature Tq = 170 MeV. This is unavoidable — for thermodynamical 
reasons the volume required for an initial temperature below Tc would be by 
far too large to be found after the expansion within a sensible thermalization 
time of 1-2 fm/c. Note that there is no mixed phase present. This is again a 
consequence of the EoS which indicates a smooth crossover between the two 
phases. Close to the critical temperature there is a soft point in the EoS and 
its effect is included in the model. 



As a test for the extracted fireball evolution scenarios, we consider dilepton 
emission. Once the thermal emission rate — ^ d^xd A p'^ >T ^ & hot source is 

known, the experimentally measured rate can be calculated as 
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(3) 



Here 'Acc' refers to the acceptance characteristic of the detector. Using 
the thermal quasiparticle model [|J for the QGP spectral function and and 
improved vector meson dominance model combined with chiral dynamics || 
for the hadronic part, we insert the resulting rate into the fireball evolution 
shown above in order to compare to the CERES data H. 
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Figure 2: Left panel: Total dilepton yield for 30% central Pb-Au collisions at 
160 AGeV (full), hadronic cocktail (dotted) and QGP contribution (dashed). 
Right panel: Contribution of the vector meson channels p, uj and 0. 

The 160 AGeV scenario is shown in Fig. 0, the 40 AGeV scenario in Fig. |3|. 
We calculate the full rate and also the contributions from different processes, 
namely the thermal rate from the QGP, the so-called cocktail contribution 
from Dalitz decays after freeze-out and the vector meson channels. In both 
scenarios we find good agreement with the data within errors. The driving 
force of the dilepton excess around 500 MeV invariant mass is the broadening 
of the p-meson due to finite baryon density. The QGP contribution is visible 
at invariant masses above 1 GeV in the 160 AGeV scenario and negligible for 
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Figure 3: Left panel: Total dilepton yield for 30% central Pb-Au collisions at 
40 AGeV (full), hadronic cocktail (dotted) and QGP contribution (dashed). 
Right panel: Contribution of the vector meson channels p, uo and <j). 

40 AGeV. Note that we do not aim for a best fit. For example, in the invariant 
mass region around 200 MeV, our calculation overshoots the data. This could 
be easily cured by raising the freeze-out temperature which corresponds to 
a reduction of the hadronic rate. However, as we choose to fix the fireball 
evolution before with hadronic observables, this is not an option any more. 

It is clearly an interesting and challenging task to also apply the fireball 
evolution model to other processes, e.g. direct photon emission or J/\l/ and see 
whether further agreement can be found, giving support to the overall scenario. 
This is work currently in progress. 
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